Fluctuation Dissipation Theorem and The Dynamical 
Renormalization Group 



Achille Giacometti'- 1 '' , Amos Maritan 1 ^, Flavio Toigo® and Jayanth R. Banavar^ 4 ) 
W Institut fur Festkorperforschung des Forschungszentrums Jiilich, 
Postfach 1913, D-52425, Jiilich, Germany 
^International School for Advanced Studies 
and Istituto Nazionale di Fisica Nucleare INFN , 
via Beirut 2-4, 34014 Trieste, Italy 
( 3 ) Dipartimento di Fisica dell'Universitd and 
Istituto Nazionale di Fisica della Materia INFM, 
via Marzolo 8, 35100 Padova, Italy 
^Department of Physics and Center for Materials Physics, 104 Davey Laboratory, 
The Pennsylvania State University, University Park, PA 16802 
(February 1, 2008) 

Abstract 

The relation between a recently introduced dynamical real space renormaliza- 
tion group and the fluctuation-dissipation theorem is discussed. An apparent 

incompatibility is pointed out and resolved. 
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Over the years, Langevin dynamics has been proven to be a successful theoretical ap- 
proach to non-equilibrium problems such as critical dynamics [|1], growth processes ||, and 
interface dynamics [Q. In a previous paper Q an exact dynamical real space renormalization 
group analysis of Langevin dynamics derivable from a Gaussian field theory was presented. 
The method has the same range of applicability as the static counterpart, namely, one- 
dimensional lattices and a whole set of hierarchical or self-similar structures. Apart from 
the methodological interest of the scheme per se, the results of the method are expected to be 
of practical interest as well, in view of the realization of growth processes on electrochemical 
electrodes |J. Unlike other real space approximate schemes @, one important feature of 
this approach is that it contains the static analog as a particular case. This naturally raises 
the question of the compatibility of the scheme with the standard fluctuation-dissipation 
theorem which relates the mobility (assumed unity in the rest of the paper) to the vari- 
ance of the noise which drives the system. In this note, we will address this issue within the 
framework of an exactly solvable model. This model, although simple, is important since it 
is the zeroth order one in any perturbative expansion. 

Consider a system described by a Gaussian field theory in 1-d with Hamiltonian (or 
action) H({<p}) given by: 

x,y 

with the field variables ip x defined on the sites x of a lattice. As it is well known ||, the 
implementation of an exact static real space renormalizaton group analysis is possible only 
in the case when the matrix A x>y has a nearest-neighbour restriction: 

^-x,y Q>x&x,y &\x— y\,l (2) 

The renormalization procedure has three basic steps : 

1) The set E of lattice points is divided into E s (sites which survive the decimation) and Ed 
(sites which are decimated) so that their union is E and their intersection is null. The 
degrees of freedom corresponding to Ed are then eliminated either by integration over 
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the corresponding field variables, or by elimination of the corresponding equations of 
motion (in the dynamical case). 

2) The surviving fields are rescaled so that the new Hamiltonian has the same form as the 

original one. 

3) Lengths are scaled in such a way that the original lattice constant is recovered. 

As result of this procedure, a mapping between the original and the new set of parameters 
is obtained. 

Let us now recall the results of the same procedure in dynamics Q. The simplest 
Langevin dynamics can be constructed as: 

where the stochastic noise rj x {t) is chosen from a Gaussian distribution which has a zero 
average and variance 

< Vx(h)vy(t 2 ) >= 2D x>y 5{t x - t 2 ). (4) 

The renormalization scheme works along the same lines as in the static case, but with the 
important difference that the renormalization of the noise has also to be taken into account. 
It was shown in ref Q that a necessary and sufficient condition for the implementation of 
the renormalization procedure was that the matrix D x>y has the form: 

D 8 Xi y + Di5| a ._j,| ) i. (5) 

Since both D and Di are different from zero under renormalization, the minimum param- 
eter space for the fluctuation matrix D is given by {Dq,D\}, implying nearest-neighbour 
correlation of the noise. We stress once again that both the static and the dynamic schemes 
we are dealing with are exact, i.e. closed in the parameter space. 

We now turn to the connection with the fluctuation-dissipation theorem. The Fokker- 
Planck equation associated with the Langevin dynamics (^[) with a noise whose variance is 
(1) is 0: 



§-P(W}, t) = £ ^-{P(W},t)^-H({ V }) + £ D * P(M, f)]. (6) 

It is then easy to show from (||) that, if we denote by Peq({y?}) the equilibrium probability 
distribution obtained in the t — > oo limit, a necessary condition for (A/" _1 is a normalization 
factor) 

Pe^})=M- l e~ H ^ (7) 

to be satisfied is that the matrix D x y has a diagonal form, so that in the case of eq. ([5]) 
Z?i = 0. Indeed @, (f|) and ([5]) with Di = lead to the equilibrium distribution given 
by (0). However the dynamical renormalization group violates , in general, the fluctuation- 
dissipation theorem since one has to start with a matrix with a non-diagonal form and the 
equilibrium distribution of (7) is not obtained. 

We now argue that this inconsistency is only apparent. We will show that the D\ ^ case 
will lead to an equilibrium distribution of a more general Gaussian Hamiltonian. However, 
we will explicitly demonstrate that equilibrium correlation functions corresponding to the 
two Hamiltonians are characterized by the same spatial decay. 

Let us consider a more general Gaussian Hamiltonian: 

^(M)=nE^4»^» ( 8 ) 

Z x,y 

where the symmetrix matrix A x y is not necessarily restricted to nearest-neighbours. On 
again using eq. @, it can be shown readily that a Langevin equation of the more general 
form: 

leads to the equilibrium distribution: 

Peq(M)=AT- 1 e-^. (10) 

Thus if one starts with the model defined by eqs. (|I|-(|) with D = 1 and Dx — 0, so 
that the equilibrium distribution is given by eq. (^), then the Dynamical Renormalization 



Group (DRG) leaves the Langevin equation of the form (0) with A and D given by eqs. (0) 
and Q) respectively, with Di ^ after one RG step. 

From eqs (8-10), the equilibrium distribution is given by (|10|), with H given by (pi) and 



A = D~ X A. (11) 

A is symmetric since both D and A depend only on the difference \x — y\. Thus initially 
A = A and then the DRG flow for H occurs in a wider parameter space than the original 
one given by eq. (g) (i.e. A no longer has the form @ ). On the contrary a static RG would 
leave H of the form (HH^). However, as we show below, these two Hamiltonians are equivalent 
in the sense that the corresponding correlation functions have the same leading behaviour. 
It is noteworthy that a static Gaussian model with interactions defined by A appearing in 
the Hamitonian (j8|) not restricted to short range, can be exactly renormalized through the 
dynamics, provided that the matrix D ■ A is restricted to nearest-neighbor interactions! 

As an explicitly solvable example, let us consider the simplest case of a one dimensional 
lattice with a x = a. The static recursions are [HI: 



a' = a 2 - 2. (12) 

The corresponding dynamics leads to the following recursions Q: 

c*V) = aV) - 2 (13a) 
D , Q =^D + D 1 + o(u) (13b) 
D[ = ]-D + \d, + o(u) (13c) 

o I 

where a(uS) = a — iu and to is the Fourier variable conjugate to time. 



From equation ( |13a|) it is apparent that the static case is recovered in the limit oj — » 
(i.e. t — * oo), that is the statics is included in the dynamics as we mentioned already. 
Let us now discuss the two-point correlation function given by: 



n 

x >y 



<<Px<P y >=^J V V e -^M) Wy = (A- 1 )^ (14) 



where Z is the partition function and T>ip = Yldip x . A similar expression holds for G x<u 
corresponding to the Hamiltonian H({ip}). 

Due to the particular form of A and D and using eq.flTID, one can easily see that the 
G Xj y and G Xj y are given by (with the lattice constant set to unity): 

G X ,y= F i (15) 

J-tt 2n a-2cosg 

and 

G xy = r^e^ Do + 2DlC ° Sq . (16) 
J-tt 2ix a — 2 cos q 

The integrals can be carried out exactly and in the a > 2 case, where the correlation functions 
are real, they yield the same coarse-grained behaviour i.e: 

G -=VH^?"*"" (17) 



(when \x — y\ > 1) and 



e„ = R±±E^ e -M*-'A (18) 

\Ja l — 4 



where we have defined: 



A=mf(i-yi-^)]i (19) 

This shows that the two Hamiltonians H and H are indeed equivalent, in the sense that the 
corresponding correlation functions decay similarly. 

Conversely, one may show that A may be written in the form (|ll|) with A and D given 
by (0) and @ respectively, provided that, at leading order, it is of the form: 

4 B)V = a8 x , y + be-'4*-v\. (20) 

In summary we have discussed the issue of the relationship between the dynamical renor- 
malization group and the fluctuation-dissipation theorem, which was prompted by our pre- 
vious analysis. We have showed that although at first glance there is a violation of the 
fluctuation-dissipation theorem, a more careful analysis restores its validity. 
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